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Abstract 

Spherically symmetric solutions of the SU(5) Einstein- Yang-Mills-Higgs system 
are constructed using the harmonic map ansatz [1] . This way the problem reduces to 



solving a set of ordinary differential equations for the appropriate profile functions. 

hj ; 
m ■ 
<N : 1 Introduction 

>• ■ Magnetic monopoles are of diverse interest since they are predicted from grant unified 

theories (GUT) and embody a rich mathematical structure. Also, they appear in non- 

{Nj \ perturbative field theories and provide a new perspective on particle physics phenomenol- 

ogy. In particular, the SU(5) gauge group plays a central role in GUT and thus it is 
natural to classify the magnetic monopoles related to this model [2] . A few years ago, the 
effects of gravitation on monopoles were (also) considered [3] and revealed a rich pattern of 

,-G solutions (including the occurrence of black holes) related to the gravitational parameter: 

a 2 = Gv where G denotes Newton constant and v G K is the vacuum expectation value 
of the Higgs field. More recently, a new interest for SU(5) monopoles was stimulated by 
the discovery of a deep analogy between their magnetic charges and the electric charges in 
one generation of elementary particles [4]. This originated several new papers on the topic, 

S^ see e.g. [5, 6] and references therein. Here we use the harmonic map ansatz [1], recently 

applied to SU(3) gravitating monopole [7], in order to construct their SU(5) counterparts. 
A similar analysis has been applied for deriving SU(5) solutions (including black holes) 
which are embeddings of the SU(2) ones [8] . However, the solutions constructed here are 
non- embedded of the SU(2) ones and correspond to monopole-antimonopole configurations. 
The SU(5) Einstein- Yang-Mills-Higgs action is given by: 



S 



\tr (F^ F^) - tr (£>„$ £>"$) - U($) 



^d A x (1) 



16nG 2 
where the potential is of the form [4] : 

U($) = -Axtr^ 2 ) + A 2 (tr($ 2 )) 2 + A 3 tr($ 4 ) - V mm . (2) 

Here g denotes the determinant of the metric while the field strength tensor is defined 
bv ; F^ u = d^A u — d v A^ + [Afj,, A u ] and the covariant derivative of the Higgs field reads: 
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D^ = <9 M $ + [A^, $]. The matrix 77 represents a constant matrix of the form: 77 = wIn, 
where ljv denotes the unit matrix in N dimensions. Finally, V min = — 15Af /(60A 2 + 14A 3 ) 
has been subtracted due to the finiteness of the energy. 

The boundary conditions are such that the energy is finite and the Higgs field at infinity 
is a given constant matrix $(0, 0, 00) = i$ 

$0 = diag(«i,« 2 ,«3,«4,«5) (3) 

in a chosen direction (while since $ G su(5): Yl\ K % — 0). In addition, the asymptotic 
values Go of the magnetic charge G — (1 + \z\ 2 ) 2 F Z z is given by 

G = G(0, 0, 1) = diag (n 1 , n 2 - n 1 , n 3 - n 2 , n A - n 3 , -n 4 ) . (4) 

Variation of (1) with respect to the metric g^ v leads to the Einstein equations 

Ri*u - 7:9 ^uR = 8vrG T^v (5) 

with the stress-energy tensor T^ u = g^L — 2 J^ given by 

T^ = tr (2D^ D u $ - g^D a § D a $) + 2tr [g af3 F^F^ - \g^ F afS F a ^ - g^V^). (6) 

In what follows we consider the static Einstein- Yang-Mills-Higgs equations in order to 
construct their SU(5) spherically symmetric and purely magnetic (ie A = 0) solutions 
based on the harmonic map ansatz first introduced in [1]. 

2 Spherical Symmetry 

The starting point of our investigation is the introduction of the coordinates r, z,z on 1R . 
In terms of the usual spherical coordinates r, 0, the Riemann sphere variable z is given 
by z = e"Han(#/2). In this system of coordinates the Schwarzschild-like metric reads: 

4r " 2 j_ J= D /^ _ 1 2m ( r ) 
B(r) ' (1 + UI 



ds 2 = -A 2 {r)B{r)dt 2 + -^dr 2 + - ' dzdz, B{r) = 1 - tlH^LL^ (7) 



where A and B are the metric functions which are real and depend only on the radial 
coordinate r, and m(r) is the mass function. The (dimensionfull) mass of the solution is 
given by m^ = ra(oo). For this metric the square-root of the determinant takes the simple 

form: 

2r 2 



-g = iA(r) 



(1 + M 2 ) 2 ' 

Then, the action (1) simplifies to 

S =/ { *^tr Q*.l') + ^tr (F$ -BU «*.)') - (l±ffitr (| D ..|>) 

-V($)X^gr 2 drdt (9) 



and the matter equations can be obtained by its variation with respect to the matter fields. 
In addition, the Einstein equations (5) take the form: 



— m' = 8vrGT °, 



-—B = &kG (T*-T r 
r A \ » r 



(10) 



where prime denotes the derivative with respect to r, and 

B(l + \z\ 2 ) 2 



T °= (1 + I ! |2)4 trirr 
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-t 



i(\F rg \ 2 )-Btr((D r ^ 



[l + \z 



2\2 



tr(|L> z $| 2 ) -V($) 



T n u - TI 



25(1 + Id 2 ) 2 



tr(|F rz | 2 ) -2-Btr((A.$) 



Next we introduce the harmonic map ansatz for the Higgs and gauge fields [1] 
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where hj(r), gj{r) are the radial depended matter profile functions and P(z, z) are N x N 
Hermitian projectors: P 3 = PJ = P 2 , which are independent of the radius r. Note that all 
N — 1 projectors Pi are orthogonal to each other since P%Pj = for % ^ j and that we are 
working in a real gauge, since A z = —A\. As shown in [1], the projectors Pk defined as 



Pk 



(A k ffA k f 



|A fe /| 



fc -f 12 



k = 0,..,N-l 



(13) 



where A/ = d z f- 



■ | f i/ give the required set of orthogonal harmonic maps (for details see 
). Moreover, the spherically symmetric harmonic maps can be constructed by applying 
the orthogonalization procedure to the initial holomorphic vector 



/= (l,2z,V6z 2 ,2z 



3 „4\f 



(14) 



Then, under the transformation: hj = J2t=j bk and c 3 - = 1 — g 3 — g 3+ \ for j = 0, . . . , 3 and 
<74 = 0, the equations of the profile functions b ,bi,b 2 ,b 3 and c ,Ci,c 2 ,c 3 can be obtained 
from variation of (9). In fact, the energy-momentum tensor Tg can be evaluated explicitly: 
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Co 2 + \ (C 2 + c'i) + c; 
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4K + 2 ( c i & i + 4H) + 4H 



3b[ , b> 3 



b'o 2 + \ (b\ 2 + b 12 ) + b'i + b l o[b l 2 + -f + -^)+ b' 3 (b\ + ^ 



3 ,/ 



26' 1 & , 2 



c* + 18 c? + c* + 8c: 



. 3 - 4cg - 6 [cf + 4)-Aci- 12 c^ + c z 2 c 



2^2 
3 



18cfc£ + 10 



vm 



(15) 



where 



vm-f 



«+§(«+«) +n+Ut,+^+fj+h,h l +^ 



26i6 2 
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A3 



^ + ^(^ + ^)+^ + &o(fe 2 + ^ + |)+&3^ + ^)+2M 2 ] 2 



125 



52 



&o (6 3 + 36i+26 2 ) +&1 (6o+36 2 +26i 



-56 



6? {h+^+2b 



■b\ L+^W 



+42 (b\ + 6*)+52 (6* + b£) + 12 [6 6i (116 + 7&i) (63 + 26 2 )+6 2 6 3 (II63 + 76 2 ) (6 + 26i)] 
+1926 2 6i (6 2 + 63) (6 + 61) + 108 [b 2 b 2 (6 2 + 63) + 6^61 (61 + 6 )] + 198 (b\bl + 6^6?) 

+426^1 -V mm . (16) 



It can be seen that the energy is finite providing the functions approach their asymptotic 
values at least as fast as 1/r, and if (in addition) the constraints: Cj (00)67(00) = are 
imposed, for all j. 

In order to read off the properties of a given solution we need to compute the Higgs 
field and magnetic charge at z — 0. Explicitly, these are given by 



(17) 



$0 = - diag[6 3 +36i+26 2 +46o, 6 3 + 36 1 +26 2 -6 , 6 3 -6 -26 1 +26 2 , 
63 - 26i - 36 2 - 6 , -46 3 - 26i - 36 2 - 6 C 



Go = diag 4 1 - d , 2 1 + 2d - 3d ± ) , 6 [c( 



2 > 



-2(1 + 2d -3d), -4 1- 



which (also) determine the boundary conditions of the matter profile functions. 
After some algebra, it can be shown that the Higgs profile functions satisfy the 
ordinary differential equations: 



4)) (18) 
following 



1 
A 

1 
A 

1 
A 

1 
A 



(ABc> o y 



bo c o + -r Co (4cq - 3c? - 1 



/ v 



(AB c[ 



(ABc> 2 )' 



bic 1 + —c 1 (6d 1 -2d -3d 2 -l 



^2 c 2 + -5- c 2 (6c 2 - 2c 3 - 3c x - 1 
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^3 c 3 + -^ c 3 (4c^ - 3c: 



while the profile functions of the the gauge fields satisfy: 



{r 2 ABb' Q )' 
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bl + l(bl + bl) + bl + b (b 2 + ^ + fj + b 3 {b 1 + ^j+2b 1 b 2 
136^ + 27bj + \2b\ + 3b 2 3 + 186 [h - -^ + ^ + I861 (63 + 26 2 ) + 126 2 6 3 



(^ABb^)' 1 



2Ar 2 



— 66iq - 36 2 c 2 - 26 c - -A1&1 



1 



4Aj 
IT 

A3 

25 



■61 



36o 



36/ 



2&!& 2 



61 



&G + 2^ + ^J+^ + &3^1 + ^ + yJ+&o(fc + , ( 
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(66 2 ^ - 3b lC j - 2b 3 cfj - ^X,b 2 

68 + |(6? + ^)+^ + 6b(t2 + ^ + |)+6s(6i + ^) +26162 
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25 



36" + 126? + 76^ + 36^ + 66 26 x + -j - 6 3 - 36 3 (4&i + 6 2 ) + 661 6 2 
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Finally, the Einstein equations (10) take the form: 
— m' = 8vrGT °, 



I A' f 4 

= SvrG - 



/2 1 u / /2 1 /2\ 1 /2 
C + 2 ( C l + C 2 J + C 3 
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6o 2 + ^(6 , i 2 + 62 2 )+63 2 + 6 (6 2 



/ (u , 3 A + ^ 
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(21) 



6' 3 ( 6; + ^6' 2 ) + 2b\b> 2 



f22) 



(20) 



where ?7i(r) and Tq are given by (7) and (15-16), respectively. 

The above system of equations has to be solved with specific boundary conditions 
which ensure the regularity of the solutions and the finiteness of the ADM mass defined as: 
Madm = m(oo)/a 2 . The Einstein equations imposes the following boundary conditions 



for the metric functions: m(0) = and A(oo) = 1. The latter condition fixes the invariance 
of the equations under the arbitrary scale A(r) — > k A(r) (for k constant) and implies that 
space-time is asymptotically flat. On the other hand, the regularity of the matter fields at 
the origin requires Cj(0) = 1 and bj(0) = while the finiteness of the ADM mass implies 
bj(oo)cj(oo) = 0. However, the specific choice of the boundary conditions on bj(oo) and 
Cj(oo) is determined by the type of solution (e.g. maximal or minimal symmetry breaking) 
we are interested in. 

It is worth mentioning that, in absence of potential, the "length" of the Higgs fields is 
not fixed since when $ — > A$ and r — > r/X the ADM mass scales according to 

M ADM {\$) = X M ADM {$). (23) 

This is true also in the flat limit (i.e. for a = 0) where the ADM mass is interpreted as 
the classical energy of the solution. 

3 Numerical Results 

3.1 Maximal Symmetry Breaking Solutions 

First, we discuss solutions with maximal SU(5) symmetry breaking that is when all the 
eigenvalues of $o (or any permutation) are different. Since there are many possibilities (in 
fact, 120 possible permutations exists) we limit ourselves to few generic cases. 

The simplest case corresponds to the self-dual (SD) solution (i.e. solution of the Bo- 
gomolny equations) where $o — diag(2, 1, 0, — 1, — 2) implying that the monopole masses 
are equal to unity since 6o(oo) = &i(oo) = 62(00) = 63(00) = 1 while the gauge functions 
Cj(r) vanish asymptotically. Then, the magnetic charge is Go — diag(4, 2, 0, — 2, — 4) i.e. 
(ni, n 2 , 713, 714) = (4, 6, 6, 4) and the corresponding mass is equal to M^ x = X^=i bjUj = 20. 

Another choice would be: 60(00) = 63(00) = 3, 61(00) = 6 2 (oo) = —2, which corre- 
sponds to the non self-dual solution (NSD) where G = diag(2, —4, 0, 4, —2) implying that 
(ni,n 2 , n 3 ,n 4 ) = (2,-2,-2,2). As expected, the solution cannot be constructed analyt- 
ically since the corresponding equations are not integrable. However, it can be obtained 
numerically and its mass is evaluated to be equal to M^^ = 27. 

When these solutions are coupled to gravity (i.e. «^0), numerical simulations indicate 
that they deformed to "gravitating" SU(5) monopoles while their presence progressively 
deforms space-time. For instance, the function B(r) develops a minimum B = B m at some 
intermediate value of r = r^, as a increases. At the same time, A(r) takes its minimum 
value at the origin while at infinity tends to the value A(oo) = 1. The metric functions A(0) 
and B m decrease as a increases for the self-dual (line SD, MAX) and non self-dual (line 
NSD, MAX) solutions, as indicated in Figure 1. Similarly, in Figure 2 the a dependence of 
the product olMadm is plotted (using the same conventions for the various lines) and shows 
that the ADM mass and the product olMadm decreases and increases (respectively) as a 
increases. It should be stressed out that, due to the peculiar normalisation of the Higgs 
field, the energies are not directly comparable; therefore, the ratio o:Madm/\&o\ should be 
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considered instead. Both branches stop at some maximal value of a ; i.e. the self-dual 
solution can be deformed by gravity up to a m ~ 0.63; while the non self-dual one exists 
up to a m ~ 0.45. As in the £77(3) case [7], the region of a in order gravitating monopoles 
to exist decreases as the mass of the flat solution increases. However this is not the end of 
the story. The main branches of gravitating solutions (i.e. the non-gravtitating ones) are 
completed by secondary branches which exist on a rather small interval of a, as seen in 
the SU(2) [3] and SU(3) [7] models. Indeed, the secondary branch exists in the following 
intervals: 

SD, MAX : a cr « 0.621, a m « 0.627 
NSD, MAX : a cr « 0.424, a m « 0.448 (24) 

and as a — * a cr the minimum of B(r) tends to zero which means (in this limit) the solution 
develops a horizon (shown in Figure 1). In fact, B m — > faster than 71(0) in terms of the 
critical value of a which means that the £77(5) gravitating monopole bifurcates into an 
extremal Reissner-Nordstrom black hole. This configuration corresponds to solutions of 
the abelian Einstein-Maxwell equations and can be embedded into the non-abelian ones. 
Its mass is equal to 

m RN = m^RN — (25) 

2r 

where Q is the charge of the black hole and can be read off from the energy-momentum 
tensor, i.e. 

^=[8(4 + 4) + is (ci + 4) - <4 + 4) - e (4 + 4) - 12 [44 + 44) - 1844 + 10 

_ (26) 

Both our solutions have charge equal to Q — v^O in consistence with the numerical 
simulations (see Figure 2). 

3.2 Minimal Symmetry Breaking Solutions 

The £77(5) symmetry can be broken into many minimal breaking patterns producing so- 
lutions with non-abelian stability group. In what follows, we present two types of such 
solutions which are invariant under the £77(3) x SU{2)/U{1) group. 

First, we investigate the self-dual solution where $o = diag(3, 3, —2, —2, —2) i.e. for 
6o(oo) = b 2 {oo) = 63(00) = and 61(00) = 5 while the Cj fields satisfy the following 
asymptotic values cq — \, c\ — 0, C2 = 4g and C3 = 4=. The corresponding solution has 
energy equal to E = 30 with mass (or classical energy in the flat limit) lower compare to 
all types of solutions we investigated - when normalized appropriately. When it is coupled 
with gravity it exists up to a m ~ 0.38. Contrary to the other cases, this solution (on 
the main branch) bifurcates into a Reissner-Nordstrom solution with charge Q = \^TE as 
confirmed by our numerical analysis which does not indicate the existence of a second 
branch. Figure 3 illustrates the way the matter functions approach their (constant) values 



outside the horizon (at r = a cr ) of the approached Reissner-Nordstrom solution when 
a — 0.1 (i.e. close to the flat limit) and a = 0.3797 (i.e. close to the critical limit). 

Finally, another non self-dual solution with the same unbroken group can be constructed 
when $ is of the form $ = diag(2, — 3, 2, — 3, 2) i.e. for 6 (oo) = 62(00) = 1, 61(00) = 
63(00) = —1 and Cj(oo) = 0. The corresponding mass of the configuration is equal to 
M = 48. 

Once more, the aforementioned solutions can be considered in the presence of gravitat- 
ing fields and our numerical routines indicate that their gravitating analogues exist up to 
a maximal value of the coupling constant equal to a m ~ 0.249. In addition, a secondary 
branch exists which terminates at a ph 0.229 into an extremal Reissner-Nordstrom black 
hole of charge Q = v^O. The corresponding results are presented in Figure 1 and Figure 
2 (line SD, MIN) and (line NSD, MIN). 

4 Conclusions 

In this paper, four types of SU(5) monopoles have been constructed which can be deformed 
by gravity forming branches of solutions labelled by the gravitational coupling constant a. 
Numerical investigation of the 577(5) Einstein- Yang-Mills-Higgs equations reveals that for 
each branch, a second one exists on the interval a G [a cr , a m ] (depending on its type) in 
consistence with the results obtained in smaller gauge groups like SU(2) and 577(3). In 
fact, the solution on the second branch has a higher mass than the one with the same a 
on the first (or main) branch; while, in the limit a — > a cr , the minimum of the function 
B{r) becomes deeper and deeper and approaches zero at some intermediate value rv Ac- 
cordingly, the regular solution does not exist for a = a cr and the metric fields approach 
that of an extremal Reissner-Nordstrom black hole on the interval [r^, 00] while all matter 
fields tend to their asymptotic values. 
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Figure 1: The metric functions A(0) and B m = mm[B(r)] in terms of a. 
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Figure 2: The ADM mass in terms of the coupling constant a. 
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Figure 3: The matter profile functions of the self-dual minimal symmetry breaking solution 
for two different values of a. 
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